QCD evolution of naive-time-reversal-odd parton distribution functions 
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We reexamine the derivation of the leading order QCD evolution equations of twist-3 quark-gluon 
correlation functions, T v> f(x,x) and T^p(x, x), which are the first transverse-momentum-moment 
of the naive-time-reversal-odd parton distribution functions - the Sivers and Boer-Mulders function, 
respectively. The evolution equations were derived by several groups with apparent differences. We 

CNj ' identify the sources that are responsible for the differences, and are able to reconcile the results from 

t-H , various groups. 

o ■ 

' PACS numbers: 12.38. Bx, 13.88.+e, 12.39.-x, 12.39.St 

>V 

. Transverse spin physics has attracted tremendous attention from both experimental and theoretical communities in 
recent years [l|. These effects typically manifest themselves in various azimuthal asymmetries. The well-known exam- 
ples are the single transverse spin asymmetry in polarized proton-proton collisions [2j , Sivers and Collins asymmetries 
in the semi- inclusive hadron production in deep inelastic scattering Q, the Boer-Mulders effect in Drell-Yan produc- 
tion [4|, as well as the large cos(20) anomalous azimuthal asymmetry in back-to-back dihadron production in e + e~ 
• annihilation || . It was soon realized that these non-trivial azimuthal asymmetries in high energy collisions should be 
i^l directly connected to the transverse motion of partons inside the parent hadron. Experimental measurements of the 
I- asymmetries and the investigation to understand the underlying dynamics have provided and will continue to provide 
us new opportunities to explore QCD and the hadron structure far beyond what we have been able to achieve. 

Two complementary QCD-based approaches have been proposed to analyze the physics behind the measured 
asymmetries: the transverse momentum dependent (TMD) factorization approach [6|-|9( and the collinear twist-three 
t-H | factorization approach (Iol - [l3 |. In the TMD factorization approach, the asymmetry was attributed to the spin and 
^ ■ transverse momentum correlation between the identified hadron and the active parton, which are represented by the 
ON ■ TMD parton distribution or fragmentation function. On the other hand, in the collinear factorization approach, all 
active partons' transverse momenta are integrated into the collinear distributions, and the explicit correlation between 
the spin and the transverse momentum in the TMD approach is now included into the high twist collinear parton 
• , distributions or fragmentation functions. The asymmetry in the collinear factorization approach is represented by 
twist-3 collinear parton distributions or fragmentation functions, which have no probability interpretation, and could 
be interpreted as the quantum interference between a collinear active quark (or gluon) state in the scattering amplitude 
and a collinear quark (gluon)-gluon composite state in its complex conjugate amplitude. The relevant TMDs and the 
• • ' quark-gluon correlation functions are closely related to each other. For example, the first fci-moment of the two well- 
known naive-time-reversal-odd TMDs, the Sivers function f^ r {x^k 2 A _) [lj| and the Boer-Mulders function h^(x,k^_) 
Q, are equal to the twist-three quark-gluon correlation functions T q ^p(x,x) and T^pix, x), respectively. These two 
correlation functions are defined as [ll| 

T q , F {x,x) = J ^±^ e ™ p ^(P, ST \M0) 7 + [e STann F+{y^)] UVi)\P^t) , (1) 

T$(x >X )= f ^^ e ^rl^ (P;ST |^ (0 )[ a «+F Q +(^)]^( y r)|P, S T), (2) 

where the gauge links between field operators are suppressed and e 0123 = 1 is used. The TMD factorization approach 
is more suitable for studying scattering processes with two very different momentum transfers, Q\ ^> Qi > Aqcd, 
where the Q2 is sensitive to the active parton's transverse momentum, while the collinear factorization approach 
is more relevant for studying scattering cross sections with all observed momentum transfers hard and comparable: 
Qi ~ Q 3> Aqcd ■ Although the two approaches each have their own kinematic domain of validity, they are consistent 
with each other in the regime where they both apply [14| . 
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Both factorization approaches necessarily introduce a factorization scale, /i 3> Aqcd, to separate the calculable 
short-distance perturbative dynamics from the long-distance nonperturbative physics of the observed cross sections or 
the asymmetries. Since the physical observables, the cross sections or the asymmetries, are independent of the choice 
of the factorization scale, the scale dependence of the nonperturbative distributions [lB - l2~i| , either TMD distributions 
or twist-3 collinear distributions, must match the scale dependence of corresponding perturbative hard parts. That is, 
the factorization scale dependence of the nonperturbative distributions is perturbatively calculable and is a prediction 
of QCD perturbation theory when /i 3> Aqcd- For example, the scale dependence of the leading power parton 
distributions obeys DGLAP evolution equations whose evolution kernels are perturbatively calculable, and has been 
very successfully tested when the scale varies from a few GeV to the hundreds of GeV. 

Recently the scale dependence of the twist-three quark-gluon correlation functions T q< p{x, x) and T^j,(x, x) have 
also been studied by several groups [l8l - l2lj ]. However, there are discrepancies between these results, particularly for the 
evolution of T qyF (x, x) (also often refer to as Efremov-Teryaev-Qiu-Sterman (ETQS) function): the result in Ref. [2l[ 
is different from those in Refs. [l8l - [2Cl | by two extra terms. The purpose of our paper is to identify and resolve these 
discrepancies, and in addition we also derive the evolution equations for the other quark-gluon correlation function 

T${x,x). 
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FIG. 1. Feynman diagrams contribute to the leading order evolution kernel of quark-gluon correlation functions. 



The evolution or the factorization scale dependence of these twist-3 parton distributions is an immediate consequence 
of the QCD factorization formalism for physical observables. Since these twist-3 parton distributions are process 
independent, the perturbative evolution kernels are universal, although they could be derived in many different ways. 
In Ref. [l8j . we presented a derivation for the evolution equations of T qt p(x, x), as well as other twist-3 correlation 
functions relevant to single transverse spin asymmetries, from the perturbative variation of these functions. By 
calculating the leading order Feynman diagrams [18j . we obtained finite contributions from those diagrams in Fig. Q] 
and the following evolution equation, 



dT q! p(x,x, fj,) _ a s f 1 d£ 



din (i 2 



(I + z)T 9tF (£, x, /i) - (1 + z 2 )T q ^, & I*) 



Nc 
' 2 



where z = x/£ and P qq (z) is the splitting kernel for unpolarized quark distribution function given by 



P qq {z) = C F 



l + z z 3 r/1 
(l-z)+ 2 



and the quark-gluon correlation function Taq,f (x\, X2, n) is given by pi 



TAq,F(xi,X 2 ) 



dyi dy~ 2 
An 



ciP 2/i pi{x 2 -x 1 )P+y 2 



<P, ST |^(0) 7 + 7 5 [»4f a + (!/2)] MvT)\P,st) 



(3) 



(4) 



(5) 



The results derived in Refs. p^j2C| are consistent with ours. But, the evolution equation derived later by Braun, 
Manashov, and Pirnay in Ref. 2l| is slightly different, 



dT qtF (x,x,n) 
91n /i 2 



" (1 + z)T q AL n) - (1 + z 2 )T q , F (L £, A 1 ) 



" 2 



l-z 



- T Aq , F (x,£,(i) 



-N c S(l - z) T q , F (x, x, fx) + — [(I - 2z)T Qt p(x, x-£,(j,)- T Aq , F {x, x-£,(j,)] L (6) 
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Comparing Eqs. ([3]) and ([6]), it is clear that two results differ by two contributions listed in the third line in Eq. ([6]). 
In addition, there is a sign difference in front of the T& q F distribution in the second line. This sign difference is 
due to a fact that two groups used a different sign convention for anti-symmetric tensor g^""^: we chose e 0123 = 1. 
while Braun-Manashov-Pirnay used £0123 = 1 implying e 0123 = — ][]]. We also noticed that Ref. 22\ used the same 
convention as that in our paper, thus they obtain the same sign for the Ta 9i f term. 
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FIG. 2. Feynman diagrams contribute to the evolution from the interference of a gluon and a quark-antiquark state (left), and 
the usual interference of a quark and a quark-gluon state (right). 



In the third line in Eq. ©, the second term oc 1/2N C comes from the two Feynman diagrams in Fig. [21 The left 
diagram in Fig. [2] corresponds to the interference between a gluon and a quark-antiquark pair. This diagram was 
missed in our original calculation in Ref. ;18;j. Once we calculate this diagram, we obtain half of the 1/2N C term in 
Eq. @, 



dT q>F [x, x, n) 



<91n /j 2 



Fig. H;icft) 



2ir 



1 d£ 1 
J2N~ C 



[(1 - 2z)T q , F (x, x-£,im)+ TA q ,F(x, x - £, //)] , 



(7) 



where again the T Aq<F term has an overall sign difference due to our convention for £^ ua P . The right diagram in Fig. [2] 
is actually Fig. 7(h) in our original paper Ref. [l8j]. This diagram vanishes if the quark on the left of the cut has a 
positive momentum p + = £P + > 0, which was assumed in the original paper [l8j . However, the £ does not have to 
be larger than as long as £ + £2 > 0. By calculating the contribution from the region where £ < 0, we find that it 
gives exactly the other half of the 1/2N C term in Eq. (J6j> , 



8T q , F (x, x, fjt) 



din y? 



Fig. Hright) 



2tt 



1 d£ 1 

J2N C 



In other words, adding two diagrams together we have 



dT QyF (x, x, n) 



Fig. Ulcft+right) 2n Jx £ 2iV c 



[(1 - 2z)T q ^ F (x - £, x, n) - T Aq<F (x - f , x, (J,)] 



[(1 - 2z)T q , F {x, x-€,fi) + T Aq . F (x, x - f , //)] 



(8) 



(9) 



The other term in the third line in Eq. © , —N c T q ^ F (x,x), was missed in our original paper [18[ . The error was caused 
by a subtlety in taking the limit xi —> when we evaluate the integration J dx25{x2)x2F(x2) = lim X2 ^o ^2^(2^2) 
to get the gluonic-pole matrix element. The limit, lim^^o X2F(x2), would vanish if the function F{x2) is finite as 
X2 — > 0, which is unfortunately not always true in our calculation. We find that Fig. [IJb) and (c) have additional 
contributions to the evolution as, 



dT q , F (x,x + x 2 , /x) 



<91n /i 2 

dT qtF (x, x-\-X2,n) 



1-x 



X2 

e 2 



= it / d4 2 Tg tF (x,x + £ 2 ,m) 

Fig. [Hb) -additional * J x 2 

= ~rf \ ^2 T q>F (x + X 2 - 6, X + X2, H) 

Fig. [TJc) -additional 1 J x+xa-1 



X2 
?2 



(10) 

(11) 



Here the "■ • •" includes any regular terms which vanish safely when we take X2 — > 0. The subtlety is caused by the 
fact that the integration over d£ 2 m Eqs. (jTUJ) and (jlip is singular as X2 — > 0. To evaluate the d£ 2 integration, we first 
expand T q F {x, x + £ 2 > A*) around £2 = x 2 in Eq. (JTT 



d 

T q , F (x,x + £ 2 ,A*) = T qtF (x,x + x 2 ,y) + —T q . F (x,x + y.)k 2 ^x 2 (£,2 - x 2 ) + 

C?2 



(12) 



1 This was also speculated in Ref. [2ll . 
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The integration in Eq. (fT0|) with the first term of the expansion in Eq. (fT2|) gives 



-y / d£, 2 T q:F (x,x + x 2 ,fi) 



X2 



N Crp 

— T qtF {X,X- 



x 2 ,n) 



£2 
6 



— -Tq iF (x,X + X 2 ,fJ,) 



x 2 



,(13) 



which goes to — ifT q ^p{x, x, fi) at the limit x 2 — ¥ 0. If one assumes T q>F (x, x + £ 2 , fi) is a smooth (regular) function, 
we find that the higher order terms in the expansion in Eq. (|12[) do not contribute to the evolution in Eq. (|10[) 
when x 2 — > 0. So Fig. [IJb) gives us an additional contribution — ^ L T qtF (x, x, Similarly we find exactly the same 
contribution from Eq. Adding them together, we have 

dTq iF (x,x, n) 



... , =-N e T q , F (x,x,ii), (14) 

Crn/i Fig. [TJb+c) -additional 

which is exactly what was missed in our original paper [l8| . 

We now have a complete agreement with the Braun-Manashov-Pirnay result. In other words, in our e 0123 
convention, we have the flavor nonsinglet evolution equation for T q ^ F (x, x, fi) as 



1 



dT qiF (x,x,fi) 



2^ 



Pqq{z)Tq tF {£,£,,ll) 

1 + z 2 



N, 
' 2 



1 



(Tq, F {£,x,[j,) -T qtF (£,£,ii)) +zT,, F ((,j;,/i) +T AqtF (x,£,[j,) 



-N c 5(l - z) T qtF (x, x, fi) + [(1 - 2z)T QtF (x, x 



+T Aq , F {x,x - M15) 



Similarly, our results for flavor singlet evolution are also now consistent with the Braun-Manashov-Pirnay result. 

Using the same techniques, we could also derive the evolution equation for the other twist-three quark-gluon 
correlation function T q p(x, x, [i). The calculation is straightforward, and the result is 



9T^ F {x,x,fi) 
dhi/u, 2 



2?r L el 



A T P qq (z)T^,t^) + ^- 

N c 5(l~z)T ( q %x,x^) + 



2T^x^)-2zT^ F (^^) 



1 

2Nr. 



2(l-z)T$(x,x 



Mi 



q,F\ 



where ArF m (z) is the splitting kernel for the quark transversity given by 

A T P qq (z) = C F 



J^z)l + 2 5{l - Z \ 



(16) 



(17) 



This evolution equation was first derived in Ref. [19], which contains only the first line in Eq. (TTB1) . The first term 
in the second line, —N c T q p{x,x,iS), has exactly the same origin as those in Eqs. (fTU|) and dTTT) from calculating 
diagrams in Figs.QJb) and (c) with a caution of taking the limit x 2 — > 0. The second term, oc 1/2N C , is again due to 
the fact that the Feynman diagrams in Fig. [2] were not included in the calculation of Ref. [19| • 

In summary, we have rederived the evolution equations for both T qiF (x, x, (i) and T q J,{x, x, (i). We resolved the 
discrepancies in the literature for the evolution of ETQS function T qtF (x, x, /i). We understand that such discrepancies 
were also resolved by the other two groups 0, [24| through careful reexaminations of their original derivations in 
Refs. [ill and [2(|, also in Ref. [25| from a different approach. Using the same techniques developed in the current 

paper, we updated the calculation for the evolution of T q a F {x, x, /i) and found two additional contributions which are 
missing in the literature. These results will have important consequences, e.g., in the study of QCD resummation for 
the spin- dependent observables [l6j . 
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